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The temperature dependence of the electron spin g factor in GaAs is investigated experimentally 
and theoretically. Experimentally, the g factor was measured using time-resolved Faraday rotation 
due to Larmor precession of electron spins in the temperature range between 4.5 K and 190 K. The 
experiment shows an almost linear increase of the g value with the temperature. This result is in 
good agreement with other measurements based on photoluminescence quantum beats and time- 
resolved Kerr rotation up to room temperature. The experimental data are described theoretically 
taking into account a diminishing fundamental energy gap in GaAs due to lattice thermal dilatation 
and nonparabolicity of the conduction band calculated using a five-level k-p model. At higher 
temperatures electrons populate higher Landau levels and the average g factor is obtained from a 
summation over many levels. A very good description of the experimental data is obtained indicating 
that the observed increase of the spin g factor with the temperature is predominantly due to band's 
nonparabolicity. 

PACS numbers: 72.25.Fe 71.70.Ej 72.25.Rb 



I. INTRODUCTION 

The temperature dependence of the spin g value of 
electrons in GaAs has been a subject of controversy since 
1995, when it was shown that the experimental data, 
exhibiting an increase of the value of g as a function of 
temperature, are in a qualitative disagreement with the 
k-p theory, if one takes into account an experimental 
change of the fundamental energy gap [1]. The seeming 
qualitative disagreement between the experiment and the 
theory was confirmed by subsequent publications [2, 3]. 
The temperature dependence of the spin properties of 
electrons in GaAs is not only of academic interest since 
for possible spintronic applications the behavior of g near 
room temperature is clearly of great importance. It was 
recently shown that one can reach at least a qualitatively 



* e-mail address: zawad@ifpan.edu.pl 

** Present address: University of Konstanz and Center for Applied 
Photonics, D-78457 Konstanz, Germany 



correct description of the experimental g values if one 
includes in the k-p theory not the complete temperature 
change of the fundamental energy gap but the change 
due to lattice dilatation alone [4]. A similar result was 
obtained for the g value of electrons in InSb indicating 
that the better description was not fortuitous. 

The experiments [1-3] and their analysis presented in 
Ref. [4] have two shortcomings. On the experimental 
side, the g value data was obtained by time-resolved pho- 
toluminescence alone, while it would be desirable to ver- 
ify these data by some other method. On the theoretical 
side, the analysis in Ref. [4] is somewhat tentative as 
it does not account for the fact that at higher tempera- 
tures electrons populate many Landau levels, so that the 
measured spin g value represents an average over many 
transitions. Recognizing the importance of the subject 
and its controversial character the present work tries to 
overcome the above shortcomings in both experimental 
and theoretical aspects. First, we present new data on 
the temperature dependence of the electron g value in 
GaAs obtained with the time-resolved Faraday rotation 
and compare them with the photoluminescence data of 
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Refs. [1-3]. Second, wc complete the theoretical de- 
scription by including the electron statistics which results 
in the necessity of summation over many Landau levels 
(LLs). We show that the complete theory considerably 
improves the description of experimental data. 

The paper is organized as follows. We first describe 
new experiments on the temperature dependence of the 
electron g value in GaAs. Next, we describe the theo- 
retical procedure used to evaluate the temperature de- 
pendence of the g value and compare the theory with all 
existing experimental data. We discuss our results and 
conclude the paper by a summary. 



II. EXPERIMENT AND THEORY 

The spin g factor of electrons in a nominally undoped 
30 /im thick layer of GaAs was determined using an ul- 
trafast degenerate setup for measuring the time-resolved 
Faraday rotation [5-7]. Spin polarized electrons were op- 
tically excited by a circularly polarized pump pulse which 
was spectrally tuned to the band gap of GaAs at all tem- 
peratures. After the excitation the electron spins begin to 
precess in a magnetic field applied parallel to the sample 
surface in the Voigt geometry. By measuring the Fara- 
day rotation of a linearly polarized probe beam trans- 
mitted through the sample one is able to monitor the 
electron spin dynamics. The Larmor spin precession re- 
sults in an oscillating signal that decays exponentially 
at a rate l/Tj*, where Tj* is the ensemble spin dephas- 
ing time (Fig. 1). The g factor is determined using the 
equation of Larmor precession: g* = Iiujl/ fisB, where 
is the measured Larmor precession frequency, /Ub is 
the Bohr magnetron, and B is the applied magnetic field. 
The ultrafast femtosecond laser pulses (AA = 15.5 nm, 
Ppurnp — 10 mW, Pprobe = ImW) wcrc gcncratcd by a 
mode-locked Ti:sapphire laser oscillator operating at 76 
MHz repetition frequency. The samples were held at low 
temperatures using an optical cryostat and a magnetic 
field was gencratcid by a split-coil superconducting mag- 
net. Magnetic fields up to 6 T have been applied. 

GaAs is a medium-gap material, so that a three-level 
k-p description, successfully used for narrow gap semi- 
conductors [8, 9], is not adequate for treating its band 
structure. The reason is that in GaAs the fundamental 
gap i?o between the Fg and Fg levels is about 1.5 eV, i.e. 
it is not much smaller than the gap Ei between the Fg 
level and the upper F^ conduction level (which is about 3 
eV). It has been convincingly demonstrated that an ade- 
quate way to treat the conduction band of GaAs is to use 
a five-level model (5LM, which is equivalent to 14 bands 
including spin) in the k-p description (see Refs. [10, 11] 
and the references therein). In particular, for the descrip- 
tion of electron spin g value in GaAs-based heterostruc- 
tures the usefulness of the 5LM was demonstrated quite 
recently [12]. Since the 5LM for electrons in the presence 
of a magnetic field and its use for magnetooptical prop- 
erties of GaAs was described in some details before [10, 
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FIG. 1: Time-resolved Faraday rotation signal of nominally 
undoped bulk GaAs, recorded at a magnetic field of B = 3 T. 



11], we only mention here the main elements of this ap- 
proach. Thus the model includes exactly the Fy, Fg, Fg, 
Fy and Fg levels at the center of the Brillouin zone and 
the resulting k-p matrix has dimensions 14 x 14. There 
exist three nonvanishing interband matrix elements of 
momentum: Pq, Pi and Q, and one interband element of 
the spin-orbit interaction A. If one takes Q = and 
— (where kz is the momentum along the magnetic field 
direction) the 14 x 14 initial matrix factorizes into two 7 
X 7 matrices for the spin-up and spin-down states. These 
matrices arc soluble with the envelope functions in the 
form of harmonic oscillator functions and the eigenenergy 
problem for different Landau levels (LLs) n reduces to a 
diagonalization of 7 x 7 determinants. However, if the Q 
element is included (it comes from an inversion asymme- 
try of the GaAs crystal) the initial 14 x 14 matrix does 
not factorize and it is not soluble in terms of a single 
column of harmonic oscillator functions. Physically, this 
means that the resiilting energy bands are not spherical. 
Since the nonsphericity of the conduction band in GaAs 
is small, one can solve for the eigenenergies looking for 
the envelope functions in terms of sums of harmonic oscil- 
lator functions. This leads to number determinants com- 
posed of the fundamental 7x7 blocks on the diagonal 
coupled by nondiagonal parts involving the Q elements. 
The eigenenergies are computed truncating the resulting 
big determinants. In our computations we used typic;ally 
21 x 21 determinants. All calculations were performed 
taking the magnetic field B parallel to [001] direction 
and putting fc^ = 0. 

Now wc tmn to the temperature dependence of the fun- 
damental gap Eq. It was argued a long time ago [13] and 
confirmed by the behavior of the effective mass [14, 15] 
that the temperature dependence of the band parameters 
is governed by a dilatational change of the energy gap Eq 
and not by its total (i.e. measured) change. The other 
contribution to the change of the gap is due to phonons. 
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FIG. 2: Calculated change of the fundamental energy gap 
in GaAs due to thermal dilatation versus temperature (after 
[16]). 



However, the phonon vibrations occur on a much slower 
time scale than the time the electron needs to sample 
the interband interactions determining the effective mass 
and the g factor, which arc at optical frequencies. The 
dilatational change of the gap is given by [16] 
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FIG. 3: Spin g values for consecutive Landau levels n, cal- 
culated for T = 50 K and T = 300 K, a magnetic field B 
= 4 T and the indicated electron density N. The indicated 
points correspond to LLs that give non-negligible contribu- 
tions to the average g value. Approximate electron energies 
corresponding to LLs marked on the lower abscissa are shown 
on the upper abscissa. 



^E^\T) = -3D (^^^j\,^{T')dT' , (1) 

whore D is the bulk modulus, dE/dP is the pressure- 
induced band-gap shift, and ath{T) is the linear thermal 
expansion coefficient (see also [17]). The quantities B 
and dE/dP are readily measurable, the quantity ath(T) 
was measured for GaAs by Novikova [18] and Soma et al 
[19]. Using the values of Q!t/,(T) one calculates AEf{T). 
This was done by Lourenco et al [16], we replot their 
results in Fig. 2. It should be mentioned that athiT) 
goes through a minimum near T « 50 K, this results 
in a plateau and a flat maximum of Eo'(T) seen in Fig. 
2. The dilatational change Aii^Q* of the fundamental gap 
between K and 300 K is about 23 mcV (see Fig. 2), 
while the total change AEq°* is about 93 meV (cf. Refs. 
[21, 22]). These numbers are important when one tries 
to understand why putting the total change Ai?Q°*(r) 
or the dilatational change IS.Eq\T) into the calculations 
of g value, one obtains very different results in the two 
cases. 

Finally, we consider average values of the spin g fac- 
tor measured as a function of temperature. The mea- 
surements are usually done in relatively pure samples 
having low free electron densities. The electrons are ex- 
cited across the gap into the conduction band and into 
both spin states. The spin states are equally populated, 
but the circularly polarized light produces a well de- 
fined coherence between them. The excited electrons 
quickly thermalize and are distributed among LLs ac- 



cording to the lattice temperature without losing their 
spin or phase. Then they interfere and quantum beats 
in the photoluminescence or other effects are observed 
from many LLs. According to this picture the observed 
signal represents an average over the populated LLs in 
which the electron thermal distribution over LLs deter- 
mines their contribution to the average g value. This 
explains why various experiments performed on weakly 
doped samples give basically the same results depending 
only on the temperature. 

We assume that the electron energies are 

^nk.=^n+^, (2) 

where n is the LL number, ± signs correspond to the 
two spin states, kz is the wavevector along the direction 
of B, and TOq is the effective mass at the band edge. The 
energies £^ contain the intricacies of the band structure 
mentioned above, but we first assume a simple parabolic 
dependence on kz to make the averaging tractable. The 
spin g value is defined as (in formulas we use g* symbol) 

9* = iEX-E-kJ/t^BB. (3) 

An averaging procedure involves summation over n and 
integrations over k^ and kz- A simple calculation gives 
the average value of g* in the form 

¥{T) = ^ , (4) 
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in which the summation is over the LLs, /(£", C) is the 
Fermi-Dirac distribution function, and the square roots 
come from the integrations over kz- The integrations 
begin from the lower states £^ for each n, which can be 
either f + or £~ . 

The average g value, as given by Eq. (4), is affected by 
the temperature in two opposite ways. According to the 
fivc-lcvcl model the spin g value at the conduction band 
edge is [10], 
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where Ep^ = 2moP§/h^, Ep, = 2moP?/h^, Gq = Eq+Aq 
and Gi = Ei + Ai. The spin-orbit energies Aq and Ai 
relate to (Fy, Fg) and (r^, Tg) levels, respectively, A is 
the interband matrix element of the spin orbit interaction 
between the (F^, Fg) and (F^, F§) multiplcts (seeJlO, 
22]), and C" is due to far-band contributions. For A = 
Eq. (7) reduces to the formula given first by Hermann 
and Weisbuch [23]. 

It follows from Eq. (7) that, as the temperature T 
increases and the absolute value of the fundamental gap 
Eq decreases (see Fig. 2), the spin g^ value at the band 
edge decreases. On the other hand, with increasing T the 
electrons populate higher LLs and band's nonparabolic- 
ity comes more and more into play. The latter is known 
to make the g value less negative (see Refs. [10, 24], and 
Fig. 3). Thus, as T increases, the average g* decreases or 
increases depending on the relative strength of the two 
effects. It is now clear why putting into calculations only 
the dilatational part of the gap variation makes the first 
effect weaker, i.e. it favors increase of the g values. We 
emphasize that we do not use in our calculations Eq. (7), 
it is quoted only to make clear the dependence of (/g on 
Eq. Calculating the electron energies away from the band 
edge we deal with the effects of band's nonparabolicity 
and inversion asymmetry, as explained above. We use 
the following band parameters of GaAs at T = 0: Ep^ = 
27.865eV, Ep, = 2.361eV, Eq = -1.519 eV, Go = -1.86 
eV, El = 2.969 eV, d = 3.14 eV, C = -2.3107, C = 
-0.0375. The zero of energy is chosen at the conduction 
band edge. The above values are the same as those estab- 
lished by an overall fit of various experiments on GaAs 



(see Ref. [10]) with the exception of the far-band contri- 
bution C", as discussed below. It is assumed that only 
£^0 depends on the temperature, see discussion below. 

In Fig. 3 we plot the g values for consecutive LLs at 
two temperatures, as determined from definition (3), in 
which the energies are calculated for B = 4 T using the 
procedure described above. It is seen that, for a given 
temperature, as the energy increases with the growing 
Landau number n, the g factor increases due to band's 
nonparabolicity. On the other hand, for a given LL the g 
factor decreases with increasing temperature. The num- 
ber of points for a given temperature indicates how many 
LLs are involved in the averaging, for still higher LLs the 
occupation by electrons is so small that their contribu- 
tion is negligible. It should be noted that the g factors 
shown in Fig. 3 do not saturate at high energies at the 
free electron value of +2, as predicted by simple versions 
of the three-level or five- level k-p models (see e.g. [24]). 
The reason is that our present theory includes the bulk 
inversion asymmetry (also called the Dresselhaus effect), 
which is manifested by the appearance of the matrix ele- 
ment Q. Also, the increase of g factor with the energy, as 
shown in Fig. 3, is stronger than linear, whereas the one 
shown in Fig. 8 of Ref. [10] is weaker than linear. This is 
caused by the fact that in Fig. 3 the energy increases by 
going to high LLs (at constant B), whereas in Ref. [10] 
the energy is increased by going to high magnetic fields 
(at n = 0). 



III. RESULTS AND DISCUSSION 

Fig. 4 shows the temperature dependent g values of 
electrons extracted from time-resolved Faraday rotation 
measurements performed at i? = 4 T on undoped bulk 
GaAs (full squares). In addition, data taken from the lit- 
erature and our calculations are shown. As to the exper- 
iments, the previously published data obtained with the 
use of time-resolved photoluminescence [1-3] and time- 
resolved Kerr rotation [25] are compared with our present 
results obtained by the time-resolved Faraday rotation. 
All these time-resolved experiments measure the phase 
coherence of spin states excited by circularly polarized 
light in the Voigt geometry and differ mostly in its de- 
tection. It is seen that the obtained data agree well with 
each other and all show an almost linear increase of g 
with the temperature. However, there is some discrep- 
ancy between various measured g values at T = 0, cf. 
[23, 26, 27]. Since all the g values presented in Fig. 4 
were measured at magnetic fields i? of a few Tesla, they 
should be higher than the band edge value since they 
are affected by band's nonparabolicity. In this situation 
we take the value of the far-band contributions C", as 
given above, to give our measured value of g* = -0.452 
at T = and B = 4 T. This corresponds to the band 
edge value = -0.472. Preference to a different experi- 
mental band edge g value would simply require a slightly 
different C". 
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FIG. 4: Measured and calculated spin g values of electrons 
in GaAs versus temperature. Experimental data: full squares 

- our results, crosses - Ref. [1], triangles - Ref [2], circles - 
Ref. [3], reversed triangles - Ref. [25]. Theory: dashed line 

- calculation neglecting kz - dependence of g^; solid line - 
calculation assuming interpolated kz - dependence of g*, see 
text. 

As far as the theory is concerned, we show two cal- 
culations. The dashed Une indicates an average value of 
g computed with the help of Eqs. (4), (5) and (6), in 
which the energies are calculated for B = A T us- 
ing the procedure outUned above. It is seen that, as the 
temperature increases and more LLs become populated, 
the band nonparabolicity effect is stronger than the di- 
latational decrease of the Eq gap, and the average g* (T) 
increases. Still, the theoretical increase does not quite 
follow the experimental values for T > 60 K. The reason 
is that our simplified formula (2) assumes the spin split- 
ting to be independent of kz- This means that, as we 
average over the energy, the g value "jumps" from one 
LL to the next, as shown in Fig. 3. In reality, however, 
the g value depends also on kz- As follows from the sim- 
plified three-level k-p model, the spin g value depends 
in fact on kz similarly to its dependence on the orbital 
energy hujc{n + 1/2). In other words, the decisive quan- 
tity is the total electron energy (see [8, 24, 28]). For the 
second calculation we assume that this is the case also in 
the complete five-level k-p model and we interpolate the 
g value to vary linearly with the energy (related to kz) 
between consecutive Landau levels. If this is done, the 
averaging procedure gives the g values indicated in Fig. 
4 by the solid line. It is seen that now the theory is in 
an excellent agreement with the experimental data. The 
moral of the story is that the almost linear increase of 
the average spin g value with the temperature in GaAs is 
caused by the dominating effect of band's nonparabolic- 
ity over the dilatational decrease of the fundamental gap. 
It should be emphasized that both our experimental data 
as well as our two theoretical curves approach the abso- 
lute zero of temperature at almost vanishing slope. It 



is important that the improved experiment and theory 
exhibit this property. 

Coming to the discussion of our results, one should 
keep in mind that our theoretical description, although 
representing a considerable improvement over the proce- 
dure of Ref. [4], is still approximate. The reason is that 
we made an assumption about the kz dependence of the 
g factor. This assumption is reasonable and it allows us 
to calculate the £^ energies and carry out the averaging 
procedure over LLs in a relatively simple way. Still, in 
principle one should include the kz terms from the be- 
ginning in the computation procedure, calculate the kz 
dependence of different LLs and carry out a numerical 
summation over kz values. This, however, makes the 
problem difficult for two reasons. First, if the kz terms 
are included, one must take much larger truncated deter- 
minants to obtain convergent values of £^{kz) energies. 
Second, when performing the averaging procedure one 
must integrate numerically over kz for each LL separately 
since in principle each LL has a different kz dependence. 
Our simplifying assumption allowed us to get around the 
above difficulties. 

As to the work of other authors, Oestreich and Ruble 
[1] put into the theoretical description the total tempera- 
ture change of the energy gap but neglected band's non- 
parabolicity. Such a description predicted a decrease of 
the g factor with temperature which contradicted the ex- 
perimental results. Oestreich et al [2] and Huebner et al 
[3] included in their description band's nonparabolicity 
but neglected important features of the band structure: 
the interband matrix element of momentum Q and of the 
spin orbit interaction A between the (F^, Fg) and (Fy, 
F|) multiplets. Both these quantities have an important 
influence on the g value, which we illustrate numerically. 
Using our band parameters and A ^ the band edge g 
value is g^ = -0.472, while for the same parameters and 
A = there is = -0.346 (see Eq. (7)). The element 
Q does not enter the conduction band edge quantities, 
but at T = 300 K and B = 4 T our calculated average 
value of g including Q is 'g*{T) — -0.297, while for Q 
= the calculated average is g* = -0.396. Also, Refs. 
[2], [3] did not include the Landau quantization of the 
conduction band, so that the calculated average g value 
did not depend on the magnetic field intensity. Finally, 
and this is the decisive point, it was assumed that the 
g value depends on the temperature via the total (i.e. 
measured) change of the fundamental gap, whereas one 
should take the change due to thermal lattice dilatation 
alone. When this is done, there is no need to introduce 
a temperature variation of the interband matrix element 
of momentum Pq, as was done in Ref. [3]. Consideration 
of the momentum operator and of the wave functions 
indicates that the interband momentum matrix element 
should also scale with the lattice dilatation, its square 
should change no more than 0.1 % in the temperature 
range from to 300 K, so it can be neglected [4]. 

The latest analysis by Litvincnko et al [4]. including 
the matrix elements A and Q and the band nonparabolic- 
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ity on one side and taking only the dilatation contribution 
to the change of the energy gap on the other, was able to 
account for the first time for the experimental increase 
of the spin g factor with the temperature. However, this 
analysis assumed no kz dependence of the Landau levels 
and, consistently, took the average electron energy for 
the Boltzmann statistics to be equal to kT. This cor- 
responds roughly to our dashed curve in Fig. 4, which 
assumes no dependence of g on kz- If one includes the 
kz dependence of the energy, the average nondegenerate 
electron energy is £ = (3/2)A:T which, upon using Fig. 
3 for gn{£), corresponds quite well to the solid line in 
Fig. 4. In other words, if one applies a proper averaging 
procedure to the g value, the resulting g*{T) corresponds 
quite well to averaging the electron energy £ and taking 
g* {£) . One should add that, taking the thermal change of 
gap due to the dilatation alone, Litvinenko et al [4] were 
able to account correctly for the temperature dependence 
of the electron spin splitting in InSb. All this, together 
with the published analysis of m*{T) in GaAs [14, 15], 
strongly confirms that both orbital and spin quantiza- 
tions of the electron spectrum in a magnetic field are 
governed by the change of the energy gap due to lattice 
dilatation alone. 



spin g factor in GaAs. The time-resolved Faraday rota- 
tion technique was used to show that, in agreement with 
other data, the g factor increases almost linearly with 
T up to room temperature. This increase is successfully 
described by the five-level k-p model of the band struc- 
ture for GaAs. As the temperature increases there occur 
two effects having an opposite infiuence on the g value. 
On one hand, a decrease of the fundamental energy gap 
caused by the thermal lattice dilatation leads to a de- 
crease of g. On the other hand, at higher temperatures 
electrons populate higher Landau levels which, due to 
band's nonparabolicity, leads to an increase of the aver- 
age g{T) value. A very good agreement of our theory 
with the experimental increase of g(T) indicates that the 
band nonparabolicity effect dominates. Our interpreta- 
tion confirms the validity of the complete five-level k-p 
model for the conduction band of GaAs. 
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